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The transport properties of dense stellar electron-proton plasma is studied following an exact 
relativistic formalism in presence of strong quantizing magnetic field. The variation of transport 
coefficients with magnetic fields are found to be insensitive for the field strengths < 10 17 G, beyond 
which all of them abruptly go to zero. As a consequence, the electron-proton plasma behaves like a 
"^J " superfluid insulator in presence of ultra-strong magnetic field. 
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(N) . 1. INTRODUCTION 

The study of the effect of strong magnetic field on dense stellar plasma is one of the oldest branches of physics. It 
has gotten a new life after the discovery of a few magnetars- the strange stellar objects, with unusually high surface 
magnetic fields [1-4]. These stellar objects are believed to be strongly magnetized young neutron stars. The surface 
magnetic fields are observed to be ~ 10 15 G. Then it is quite possible that the field strength at the core region may 
go up to 10 18 G. The exact source of this strong magnetic field is of course yet to be known. These objects are also 
supposed to be the possible sources of anomalous X-ray and soft gamma emissions (AXP and SGR). If the magnetic 
fields are really so strong, in particular at the core region, they must affect significantly most of the important physical 
properties of such stellar objects and the physical processes, e.g., weak and electromagnetic interactions taking place 
at the core region. Which means, the presence of strong quantizing magnetic field at the core region should modify, 
both qualitatively and quantitatively the equation of state of dense neutron star matter, and as a consequence the 
gross-properties of neutron stars [5-8], e.g., mass-radius relation, moment of inertia, rotational frequency etc. should 
also change significantly. In the case of compact neutron stars, the phase transition from neutron matter to quark 
matter at the core region, if any, will also be affected by strong quantizing magnetic field. It has been shown that a 
first order phase transition initiated by the nucleation of quark matter droplets is absolutely forbidden if the magnetic 
O !~ field ~ 10 15 G at the core region [9,10]. However a second order phase transition is allowed, provided the magnetic 
field strength < 10 20 G. This is of course too high to achieve at the core region. 

The elementary processes, in particular, the weak and the electromagnetic processes taking place at the core region 
of a neutron star are strongly affected by such ultra-strong magnetic field [11,12]. Since the cooling of neutron stars 
are mainly controlled by neutrino/anti-neutrino emissions, the presence of strong quantizing magnetic field should 
affect the thermal history of strongly magnetized neutron stars. Further, the electrical conductivity of neutron star 
matter which directly controls the evolution of neutron star magnetic field will also change significantly. 

In another kind of work, the stability of such strongly magnetized rotating objects are studied. It has been observed 
from the detailed general relativistic calculation that there are possibility of some form of geometrical deformation 
of these objects from their usual spherical shapes [13-15]. In the extreme case such objects may either become black 
strings or black disks. In the non-extreme case, however, it is also possible to detect gravity waves from these deformed 
rotating objects. 

In a recent study on microscopic model of dense neutron star matter, we have observed that if most of the electrons 
occupy the zeroth Landau level, with spin anti-parallel to the direction of magnetic field and very few of them are with 
spin along the direction of magnetic field and Landau quantum number > 0, then either such strongly magnetized 
system can not exist or such a strong magnetic field is just impossible at the core region of a neutron star [16]. 

Motivated by the problems as mentioned in preceding two sections, in this paper we shall study the effect of strong 
quantizing magnetic field on the transport coefficients of dense stellar electron-proton plasma. We shall follow an exact 
formalism [17] which is applicable for a wide range of magnetic field strengths and obtain the transport coefficients 
from the relativistic version of Boltzmann kinetic equation by linearizing the distribution function and using relaxation 
time approximation. We shall obtain the relaxation time from the rates of standard electromagnetic processes taking 
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place inside electron-proton plasma and make necessary modification in the rate calculation due to the presence of 
strong quantizing magnetic field. We have noticed that the electrical conductivity of the medium becomes extremely 
small in presence of ultra-strong magnetic field (> 10 17 G). The magnetic field at the core region of a magnetar must 
therefore decay very rapidly (time scale ~ a few mins.) and becomes moderate or low enough. As a consequence there 
will be in principle no problem on the existence of magnetars (with very low or moderate core magnetic field). The 
formalism we have developed to obtain rates of electromagnetic processes or the relaxation time is also applicable to 
evaluate neutrino cmissivity and mean free path in presence of strong quantizing magnetic field. 

In presence of strong quantizing magnetic field, since the momentum component in the transverse plane with respect 
to the external magnetic field gets quantized, whereas the component along the field direction varies continuously 
(from — oo to +00), the momentum space volume element becomes 



d 3 p dp x dp y dp z eB n 



(2tt) 2 



(2tt) 3 



v=0 



(2 - S u0 )dp z 
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(we have assumed % = c = kj, = 1) where we have chosen the gauge = (0, 0, xB m , 0), so that the constant magnetic 
field B m is along z-direction. We have considered the simplest possible picture of neutron star matter with n — p — e 
out of thermodynamic equilibrium and the neutrinos are assumed to be non-degenerate. The baryonic components 
are interacting via a — u> — p meson exchange type mean field and the electrons are assumed to be freely moving 
particles. 

In this article, we shall first calculate the transport coefficients of electron gas. Then it is very easy to obtain the 
transport coefficients for proton matter just by replacing mass, chemical potential etc. of electrons by protons and 
taking into account the proper modification in presence of a — u — p meson exchange type mean field [6]. Spinor 
solutions for electrons in presence of strong quantizing magnetic fields are then given by 
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where | and J. represent up and down spin states respectively, 
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with H u the Hcrmitc polynomial of order v and 



(p 2 +m 2 + 2e^ m ) 1 / 2 



(4) 
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the energy eigen value with v = 0,1,2, , the Landau quantum numbers. For neutron we consider the usual 

spinor solutions. Since the temperature of the system is electron chemical potential, we have not considered the 
negative energy spinor solutions, In the energy eigen value p^_ = (2ueB m ) 1 ^ 2 is the transverse component of electron 
momentum. To overcome the serious problem on the mechanical stability and hence the existence of magnetars, we 
have studied the variation of transport coefficients, in particular the electrical conductivity of electron gas / proton 
matter in presence of strong quantizing magnetic field and try to show whether the electrical conductivity which 
is solely responsible for the evolution of neutron star magnetic field, becomes sufficiently small in presence of ultra 
strong magnetic field. 

The paper is organized in the following manner. In section 2, we have developed the relativistic version of Boltzmann 
kinetic equation for fermions in presence of strong quantizing magnetic field and obtain the expressions for transport 
coefficients following de Groot [18]. We shall use the relaxation time approximation for the collision term. In section 3 
we shall evaluate the the relaxation time from the rates of standard electromagnetic processes taking place at the core 
region. We shall incorporate the necessary changes in the rate calculation due to the presence of strong quantizing 
magnetic field (or the effect of quantized Landau levels). In the last section we shall discuss importance of the results. 
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2. BOLTZMANN EQUATION 



The relativistic version of Boltzmann transport equation is given by [18] 



(6) 



where the second and the third terms are coming from electromagnetic and gravitational (general relativistic) inter- 
actions respectively and C is the collision term. Since we have considered the flat space-time geometry and noticed 
that there is no contribution of external magnetic field, only the induced electric field arising from local charge non- 
neutrality, causes an induced current in the system will contribute in the electromagnetic interaction term and the 
curvature term is neglected. 

To obtain the shear and bulk viscosity coefficients, heat conductivity and electrical conductivity of dense electron 
gas we make the relaxation time approximation, given by 



C 



(7) 



where 



f(x,p) = fW(p)(l + X (x,p)) 



(8) 



i.e., the system is assumed to be very close to the local equilibrium configuration. Here r is the relaxation time and 
/° is the (local) equilibrium distribution (Fermi distribution) function, given by. 
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We shall now follow the general technique to obtain the perturbed part x as a linear sum of driving forces. Following 
the formalism developed in the book by de Groot [18], we express the four derivative as the sum of a space like part 
and a time like part, given by = u^D + V M with the hydrodynamic velocity, D — u^d^ is the convective 
time derivative and V M = A^d^ is the gradient operator, with A^ u = g^ v — u^u v some kind of projection operator, 
gV v = diag(l, —1, —1, —1) the metric tensor. Then using the eqns.(6)-(9), the decomposition of and the equations 
of motion, given by 
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we obtain the perturbative part x, given by 
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is the driving force for bulk- viscosity, 



is the driving force for heat conduction, 
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VP 



X ^ = W - ^A^Vvif 



is the driving force for shear viscosity and E v is the driving force for electric current (E v for v = i 
components of electric field vector) . The quantity Q is given by 



(12) 
(13) 
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Q = -|A^ + (P^ f^d - 7) + (t»(1 - A (£) - (14) 

with F(T) = P(T)/n(T), P{T) and n(T) arc the equilibrium local kinetic pressure and number density respectively. 
For a non-relativistic Boltzmann gas F(T) — T, the local temperature of the system and finally h = (e + P)/n, the 
enthalpy per particle with e the local energy density and 7 = C p /C v the ratio of specific heats at constant pressure 
and constant volume respectively. 

Now from the definition, the heat flow four current is given by 



4tH J_ oc po 



= 4 0)li + 7 9 1)AI ( 15 ) 

where the first term is the equilibrium contribution, which is identically zero. Then omitting the symbol (1), we have 
the irreversible term 
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/+00 j 
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Again using the definition 

J," - X^X qv , (17) 



we have the heat conductivity coefficient 

1 1 t /<•„. x • - /- X 

Po 

Now from the definition, the energy-momentum tensor is given by 

eB 



A = l^Xr = E( 2 - M / + J ^-r(p°u a hff^{l - /(°))A^ W (18) 



T ^ = 4^ > J 2 - M / -///(^) (19) 
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v J — 00 Po 

Which can also be written as = T^ v + T^^, with the equilibrium value 

T (o)„, = ^V(2- M / + °° ^W/ (0) (p) (20) 

and the non-equilibrium part 

T (D^ = T ^ = ^^ (2 _ M f + °°^p» x ( x , p )fW(p) (21) 

which is a symmetric second rank tensor. Now we consider a model in which a flow of electron gas with cylindrical 
symmetry is assumed. Then considering fi = r, v = z, we have from the definition 



T (i)« = (22) 
dr 



Hence the shear viscosity coefficient is given by 

^=|^E( 2 -^o) / +OC % P V) 2 / (0, (1-/ (0 V (23) 
47r 1 „ 7 -00 Po 

where p r = {2veB m ) 1 / 2 , the transverse component of electron momentum. Then it is quite obvious that in presence 
of ultra strong magnetic field, for which ^ max = 0, the shear viscosity coefficient vanishes. Here ^ max is the maximum 
value of quantum number of the Landau levels occupied by electrons for a given density and temperature. For T = 0, 
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where [ ] indicates the nearest integer less than the actual value. 

To obtain an expression for bulk viscosity coefficient, we next consider the pressure tensor, given by 

IF" = A£T CTT A^ + PA* 4 " (24) 

Where the reversible part 

n^"" = (25) 

and the non-equilibrium part 

■ E( 2 - M / +0 ° ^^KP> T x(x,P)f W {v) (26) 
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4tt 2 



Hence omitting 1, we have the traceless part of pressure tensor 
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3 4tt 2 

= JfoV^u" (27) 
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Hence we have the bulk viscosity coefficient 

* = B 2 - |A roP yrQ(l - /(0))/(0) 

We shall now calculate the the electrical conductivity for electron gas. The electric four current is given by 



(28) 



Now because of local charge neutrality, the reversible part j^ ) = and the non-equilibrium part is given by 



E(2-M/ — P*f l0) x(*,p) (30) 



4tt 2 

Then using the covariant form of Ohm's law 

= CT^^ (31) 

the electrical conductivity tensor is given by 

uv e 2 eB m ^. . . f +ao d Pz u v (v a u a 
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This equation shows that due the presence of strong quantizing magnetic field the electrical conductivity becomes a 
second rank tensor- even if the space is isotropic in nature. The ^-component is given by 



£(2 - M / 1 + ^ /<°>(1 - /<°>) (33) 



The _L_L-component is given by 
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^ = E( 2 - WueB m ) J J d Pz ^ (l + ££sl) /(0)(i - /TO) (34) 



and finally the _L z-component is given by 

f-\-OC 



* ±z = ^ = B - e -^ £( 2 " M(2-B m ) 1/2 / J (l + ^) - / W) (35) 

From eqn.(34) and (35) it is quite obvious that just like the shear viscosity coefficient, both a- 11 - and a^ 7 - components 
of electrical conductivity vanish for z/ max = 0, i.e, in the ultra strong magnetic field limit. This is of course not at all 
evident for i] v , A and a zz . 

3. RATE OF ELECTROMAGNETIC PROCESSES 

Now to obtain the numerical values of all these transport coefficients, or their variations with the strength of 
magnetic field we have to know the relaxation time t, given by 



7 = 2> 



(36) 



where Wi is the rate of i th electromagnetic process and the sum is over all possible electromagnetic processes taking 
place involving the electrons. 

Therefore to obtain the relaxation time we evaluate the rates of the basic electromagnetic processes, given by 
e + e — > e + e, and e + p — > e + p. Now in the case of e — e-Scattering it is necessary to consider both direct and 
exchange processes, whereas in the case of e — p scattering only the direct term contributes. Let us first consider the 
e — e-scattering process, then the direct part is given by 

. o 
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Similarly the exchange term is given by 
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where Q is the exchanged momentum. The rate element is then given by 

I rp{d) rp(ex) , 2 

dW = lim 1 fi fi L (39) 

t^oo t 

Since the spinors are functions of ^-coordinate through I VyPy (x) (see eqns.(2)-(4)), we use the relation 



/OO 
dx'f*(x)f(x') (40) 
-oo 



and we also need the positive energy projection operator, given by 

A+ = Y,u(k,x)u(k,x') (41) 

spin 

Substituting the positive energy spinor solutions (eqns.(2)-(3)), we have 

A+ = -^-(Ak^ifj, = and z) + mA + Bk^in = y and p y =pi_) (42) 
lhi v 
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The matrices A and B are given by 
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where primes indicate the functions of x' . Since the Dirac 7 matrices are traceless and both A and B matrices are 
diagonal in nature with identical blocks we have a few interesting relations, e.g., 



Tr(^ A.A^.B^..) = TriA^.B^.W, 
Tr( 7 ^ x ^ A 1 A 2 ..B 1 B 2 ..) = Tr(A 1 A 2 ..B 1 B 2 ..)(g^g° x - g» x g v ° + g^g v °) 
Tr (product of odd no of 7s with any number of A and/or B matrices) = 



(45) 
(46) 
(47) 



etc. The beautiful form of A and B matrices allow to multiply 7-matrices in the above relations from any side or in 
any order. The other interesting aspects of A and B matrices are 

i) k 1 ^Tr(A 1 A 2 ) = (E 1 E 2 - k lz k 2z )Tr(A 1 A 2 ) 

ii) fc lAl fc 2 ^Tr(B 1 B 2 ) = £ 1± .& J .'n(B 1 fl 2 ) 
hi) k^k^T^A^) = fc lAl fc 2 ^Tr(B 1 A 2 ) = 

iv) p 1 „k 1 ^ P2v k^Tr(A 1 B 2 ) ± = (E Vl E^ - P i z k lz )p 2± .k 2± Tv(A 1 B 2 ) 

All these results are entirely new and to our knowledge not reported before in the literature. We do believe that 
these new results can have interesting applications in various studies of properties of strongly magnetized dense stellar 
matter. With all these new formulae, the transition matrix element for e e scattering direct term is given by 
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{TriAsA^TriAiA^ + Tr(A 3 Bi)Tr(A A B 2 ) + Tr^A^Tr^^) + Tr{B z B x )Tr{BiB*)} 
+(k 3 .k 2 )(k 1 .k 4 ){Tr(A 3 A 1 )Tr(A 4 A 2 ) + Tr(A 3 B 1 )Tr(B i A 2 ) +Tr{B z A 1 )Tr{A i B 2 ) +Tr{B 3 B 1 )Tr{B A B 2 )} 
-m 2 (fc 3 .fci){Tr(^3^i)Tr(A 4 A 2 ) + Tr^B^Tr^A^} - ™ 2 Xfc 4 .fc 2 ){Tr( J 4 3 A 1 )Tr(.B 4 .B 2 ) 
+Tr(A 3 A 1 )Tr(A 4 A 2 )} + 2mi{Tr(A 3 A 1 )Tr(A 4 A 2 )}]dp 2y dp 2z dp 3y dp 3z dp 4v dp 4z 

/o(P2«)(l - /o(P3*))(l - /o(P4z)) (48) 

where j Vi = (2 - 6 Vi 0). 

The exchange term is given by 
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{Tr(A!A 4 )Tr(^ 2 A 3 ) + Tr^A^Tr^Aa) + Tr^B^Tr^B^ + Tr(Bi B 4 )Tr(B 2 B 3 )} 
+ (fc 1 .fc 3 )(fc 2 .fc 4 ){Tr(A 1 A 4 )Tr(A 2 A 3 ) + TriB^Tr^Bs) + Tr(A 1 B 4 )Tr(B 2 A 3 ) + T^B^T^B^)} 

- m 2 (fc 1 .fc 4 ){Tr(A 1 A 4 )Tr(A 2 A3) + Tr(B 1 B 4 )Tr(A 2 A 3 )} 

- m\(k 2 .k 3 ){Tr(A 1 A 4 )Tr(A 2 A 3 ) + Tr(A 1 A i )Tr(B 2 B 3 )} 
+ 2m i e Tr(A 1 A 4 )Tr(A 2 A 3 )]dp 2v dp 2z dp 3y dp 3z dp 4v dp 4z 

/o(P2«)(l - /o(P3z))(l - /o(P4*)) (49) 
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and similarly the e e-scattering mixed term is given by 
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+ 2m 2 e (k 1 .k 3 ){Tr(A 1 A 2 A 3 A A ) + Tr(BiA 2 B 3 A 4 ) + 2m 2 e (k 1 .k 2 ){Tr(A 1 A 2 A 3 A 4 ) + Tr{BiB 2 A 3 A 4 )}] 
+ 2m 2 e (k 1 .k 4 ){Tr(A 1 A 2 A 3 A A ) + Tr(BiA 2 A 3 B 4 )} + 2m 2 e (k 2 .k 3 ){Tr(A 1 A 2 A 3 A 4 ) + Tr(AiB 2 B 3 B 4 )} 
+ 2m 2 e (k 3 .k A ){Tr{A l A 2 A 3 A A ) + TriA^B^)} + 

2m 2 e (k 2 .k i ){Tr(A 1 A 2 A 3 A A )Tr(A 1 B 2 A 3 B 4 )} - m\Tr(A x A 2 A 3 A A )\ 

dp 2 ydp 2z dp 3y dp 3z dp4ydp4zfo{p2z)(l - /o(P3*))(l - fo{P4z)) (50) 
and finally the rate of e p-scattering is given by 



Wcp 



/e 4 1 
dxdx'— ^3 -5{E t +E 2 - E 3 - E 4 )5(p ly +p 2y - p 3y - p4 V )S(p lz +p 2z - p 3z - p 4z ) 
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x ^^^: [(fc3 - fc4)(fci - fc2) 

{Tr(A 3J 4 1 )Tr(A 4 A 2 ) + Tr (A 3 B^ )Tr(A i B 2 ) + Tr(B 3 4i)2V(B 4 A 2 ) + Tr(B 3 B 1 )Tr(B 4 B 2 )} 
+ {k 3 .k 2 ){k l .k 4 ){Tr{A 3 A l )Tr{A 4 A 2 ) + Tr(MB{)Tr(B A A 2 ) + TriB^TriA^) 
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-m, 



^(k4.k 2 ){Tr(A 3 A 1 )Tr(B i B 2 ) + Tr(A 3 Ai)Tr(A 4 yl 2 )} + 2m 2 e m 2 p {Tr(A 3 A 1 )Tr(A 4 A 2 )}} 



dp 2y dp 2z dp 3y dp 3z dp4 V dp4 Z fo(p2 Z )(l - /o(P3«))(l - h{p4 Z )) (51) 

To evaluate the multidimensional (eight dimensional) integrals, we have used three 5-functions and finally evaluate 
numerically using multi-dimensional Monte-Carlo integration technique and obtain the rates from the eqns.(48)-(51). 
We have generated the Hermite polynomials appearing in the matrices A and B numerically. We obtain the relaxation 
time from eqn.(36) and finally evaluate the kinetic coefficients from eqns.(18),(23),(28),(33)-(35) for various values 
of magnetic fields B m , temperature T and matter density tib- In figs.(l)-(6) we have shown the variation of kinetic 
coefficients with magnetic field strength. Identical results can also be obtained for dense proton matter. 



4. DISCUSSIONS 



From the figs.(l)-(6) it is quite obvious that the kinetic coefficients are almost independent of magnetic field for 
moderate strengths (< 10 17 G), but all of them go to almost zero for magnetic field strength beyond 10 17 G. The 
first conclusion is therefore, that at ultra strong magnetic field, the matter (electron or proton matter) behaves like 
a superfluid insulator. The mechanism of superfluidity is of course completely different from conventional neutron 
matter. 

Now as we know from classical plasma physics in presence of strong magnetic field, that the charged particles can 
only travel along the lines of force, motions are almost forbidden across the field, as a consequence <j±± or a± z vanish 
for a classical plasma in presence of ultra-strong magnetic field, whereas a zz remains non-zero even at very high 
magnetic field strength. On the other hand, in the case of quantum mechanical plasma as the magnetic field become 
strong enough, the maximum value of Landau quantum number ^ max — > 0, which further means the system becomes 
effectively one dimensional in the ultra-strong magnetic field limit. Then in the collision dominated scenario, since p z 
varies continuously from —00 to +00 and p± — > 0, a zz — > 0. To elaborate this point a bit- since the current is flowing 
along the field lines only, then from the symmetric nature of p z , we have 

J ( z +) = -3 { z - ] Jz = J ( z +) + j<"> = (52) 

This is specially true for an effectively one dimensional collision dominated quantum plasma. In the case one dimen- 
sional flow (free streaming) of charged particles along the field lines, this is not true; a zz has non-zero finite value. 
Since all the three component of electrical conductivity vanish in presence of ultra-relativistic magnetic field, this 
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must affect significantly the Ohmic decay of strong magnetic field. We know that the Ohmic decay time scale is given 

by 



where I is the dimension of the system. Since a — > for extremely strong magnetic field, as a consequence, the 
field should decay quickly. We have noticed that a field of strength 10 18 G become ~ 10 15 G within a minute. This 
moderate field of course will take longer time to decay further (long decay time scale). Therefore the last conclusion 
is that stronger the magnetic field, shorter the decay time scale and as a result we do believe that the magnetic field 
strength at the core region of a magnetar can not be high enough and consequently there will be no problem on the 
mechanical stability or the existence of magnetars. 




B / B (c)(e) 

FIG. 1. Variation of shear viscosity coefficient with magnetic field B. Curve (a): T — 15MeV and ns = 3no, Curve (b): 
T = 15MeV and hb = 5no, and Curve (c): T = 30MeV and hb = 3no, 
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FIG. 2. Variation of bulk viscosity coefficient with magnetic field B. Curve (a): T — I5MeV and ub = 3no, Curve (b) 
T = I5MeV and tib = 5no, and Curve (c): T = 30MeV and ris = 3no, 




i i i i nii| 1 i i 1 1 ii| 1 I 1 1 in 1 I 1 1 in 1 i i 1 1 1 

1 10 10 2 10 3 10 4 10 5 

B/B (0(e) 

FIG. 3. Variation of heat conductivity coefficient with magnetic field B. Curve (a): T = I5MeV and ns = 3no, Curve (b) 
T — I5MeV and ub = 5no, and Curve (c): T = 30MeV and ub = 3no, 
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FIG. 4. Variation of zz component of electrical conductivity with magnetic field B. Curve (a): T = 15MeV and n B = 3no 
Curve (b): T = 15MeV and n B = 5n , and Curve (c): T = 30MeV and n B = 3n , 
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FIG. 5. Variation of pp component of electrical conductivity with magnetic field B. Curve (a): T = 15McV and ns = 3no 
Curve (b): T = 15MeV and n B = 5n , and Curve (c): T = 30MeV and n B = 3n , 
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FIG. 6. Variation of zp component of electrical conductivity with magnetic field B. Curve (a): T = 15MeV and ns = 3no, 
Curve (b): T = 15MeV and n B = 5n , and Curve (c): T = 30MeV and n B = 3n , 
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